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ABSTRACT 

We investigate the luminosity and redshift dependence of the quasar continuum by means of the 
composite spectrum using a large non-BAL radio-quiet quasar sample drawn from the Sloan Digital 
Sky Survey. Quasar continuum slopes in the UV-Opt band are measured at two different wavelength 
ranges, i.e., 0^12 (1000 ^ 2000A) and a ,^24 ( 20 00 ^ 4000A) derived from a power-law fitting. Gener¬ 
ally, the UV spectra slope becomes harder (higher a^) toward higher bolometric luminosity. On the 
other hand, when quasars are further grouped into luminosity bins, we find that both 0,^12 and ai, 2 i 
show significant anti-correlations with redshift (i.e., the quasar continuum becomes redder toward 
higher redshift). We suggest that the cosmic dust extinction is very likely the cause of this observed 
— z relation. We build a simple cosmic dust extinction model to quantify the observed reddening 
tendency and find an effective dust density nay ^ 10~^h Mpc“^ at 2 ; < 1.5. The other possibilities 
that could produce such a reddening effect have also been discussed. 

Subject headings: dust, extinction — quasars: general 


1. INTRODUCTION 

Obser vationally, quasars s how very similar UV-optical 
spectra (|Francis et al.lfT99^ , which can be characterized 
by a featu reless continuu m and a series of broad emis¬ 
sion lines (|Petersonlll997h . At wavelengths longer than 
the Lya line, from 1300 to 5000A, the continuum of 
quasars can be well described by a power-law /i, oc 
where ai, charact erizes the continuum slop e with a typi¬ 
cal value ^ — 0.5 l|Vanden Berk et'aI1l200lD . In terms of 
more negative a^, means a steeper (or redder, softer) 
continuum. 

There is a general agreement that the UV continuum 
of higher luminosity quasars is harder (ICarballo et al.l 
Il999t Telfer et al.ll200^ fDavis et al.l 120071 however, see 
also iKennefick fc BursickI l2008h . On the other hand, 
there has been a long running debate on the existence 
of the redshift dependence (evolution) of the quasar con¬ 
tinuum. Some studies suggest the s light hardening of the 
continuum toward higher redshift (jCarballo et al.l [19991 
IKennefick fc Bur sickll200 8H. while others find hardly any 
evolution (iKuhn et al.l ^Oll : iPentericci et al.l l2003ll or 
even reverse (IWrightll198i]L 

Several mechanisms could cause the redshift depen¬ 
dence of the quasar continuum, e.g., the evol ution of 
the accretion process (jKennefick fc Bursickll200^ and/or 
the evolution of the dust component. Among them, one 
of the interesting mechanism is the intergalactic (cos¬ 
mic) dust extinction, which makes the quasar spectra 
statist ically redder at higher redshifts. Indeed, IWrightI 
(|l98lh attributed the redder spectra toward higher red¬ 
shift they detected to the effect of the cosmic dust, and 
concluded that the average line of sight has cosmic ex¬ 
tinct ion_o£j4^^= RSSi 0.51 out to z=3 (however, see 
also ICheng et al.lll99lH . Compared to other attempts 
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on probing the cosmic dust extinction (e.g.. iMore et al.l 
120091 : iMenard et akll^lOal) . quasars take advantage of 
their high luminosity and large sample size; therefore, 
they can be explored to higher redshifts. 

In this Letter, we take a l arge quasar samp le from Sloan 
Digital Sky Survey (SDSS, lYork et ^l2000h and use the 
continuum slope of the quasar composite spectrum to 
probe its redshift dependence, then further explore its 
implications for the cosmic dust. 

This Letter is organized as follows. In Section 2, we 
describe the data and method of building the composite 
spectrum. In Section 3, we report the finding of a signifi¬ 
cant dependence of a^, on the redshift: redder spectra at 
higher redshifts. In Section 4, we build a simple cosmic 
dust extinction model to quantify the reddening effects. 
In Section 5, we discuss other possibilities. Finally, we 
present our conclusions in Section 6. Throughout this 
Letter, we adopt the cosmological parameters Ha = 0.7, 
Dm= 0.3, and h = 0.7. 

2. DATA 

The SDSS legacy survey provides a database of quasar 
samples that includes photometries in u,g,r,i, and z 
bands and spectra from 3800 to 9 200 A with resolu¬ 
tions of about 2 000 (lYork et al. l2000ll. In the SDSS DR7 
quasar catalog. fSchneider et al. HoiOD compiled 105,783 
quasars that are more luminous than Mi = —22.0 and 
have at least one emission line with FWHM larger than 
1000 km s“^ or have interesting/complex absorption fea¬ 
tures. This catalog provides the basic database for our 
study. 

Accompanying the D R7 quasar catalog , there are other 
value-added catalogs. iHewett fc Wildl (|2010l) recalcu¬ 
lated the redshift of quasars and reduced the systematic 
effects of the redshift to the lev el of 30 km s“^ per unit 
redshift. IWild fc Hewetd (|2010h improved the sky sub¬ 
traction in the red wavelength band. After measuring 
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the spectroscopic features, IShen et HI (l201lD calculated 
the implied physical parameters (e.g., the bolometric lu¬ 
minosity, black hole mass, and Eddington ratio) for each 
quasar. 

In t his study, we take the spectra from lWild fc Hew^ 
(l2nTnli and u se th e improved redshifts from 
iHewett fc Wildl (120101) . We exclude the quasars 
with broad a bsorption lines fro m our analysis according 
to the flags in IShen et al.l (|2011ll . The radio-loud quasars 
show significant re dder continua than radio-quiet ones 
(iLabita et al.l l2008f) and the fractio n of them change s 
with both redshift and luminosity (j.Tiang et al.l l2007ll . 
For these reasons, we also remove them. Therefore, 
our final sample with calculated bolometric luminosity 
includes 91,131 objects. 


2.1. Composite spectrum 

We hrst correct the galactic reddening f or each quasar 
using the SFD mao (ISchlegel et ahlllQQSH and the red¬ 
dening curve of ICardelli et al.l ( 1989lb Then, we mask 
bad pixels and shift the s pectra to their r est fr ame ac¬ 
cording to_thej;edshifts_of|^wt^^Wi3 ([2013) ■ Next, 
we follow IVanden Berk et al.l ( 20011) to make the quasar 
composite and calculate the statistical error for the com¬ 
posite spectrum by dividing 68% semi-quantile to the 
square of the contributing number at each wavelength. 

The composite spectrum of all 91,131 quasars is shown 
in Fig. [T] as t he blue line. For comparis on, the composite 
spectrum of IVanden Berk et al.l (|20011) is shown as the 
red line. As can be seen, these two composites show good 
consistency. Because of a much larger sample size, our 
composite spectrum has a much higher signal-to-noise 
ratio (shown in the sub-panel of Fig. [IJ. 


1350-1365 A and 4200-4230 A, have been adopted, which 
avoid the Lya line at shorter wavelengths and contamina¬ 
tion from the host galaxy at longer wavelengths. How¬ 
ever, this range is too long to be covered by any sin- 
gle SPSS s p ectrum at any given redshift. Similar to 
iDavis et"^ (l2007l) . we add a new continuum window 
2210-2230 A. With this configuration of continuum win¬ 
dows and the 3800-9200 A wavelength range of SDSS 
quasars, the far-UV slope a^-i 2 and the near-UV slope 
0,^24 can be htted for quasars in the redshift intervals 
1.80 < 2 < 3.15 (n = 22251) and 0.71 < 2 < 1.19 
{n = 18068), respectively. 

3. RESULT 

3.1. Trend of ai, with bolometric luminosity 

We bin the quasars into four bolometric luminosity 
(Lboi) bins with a bin size of 1 dex from 10^^ erg s“^ 
to 10^® erg s“^ and then make composites respectively. 
To show the possible variance of the composites more 
clearly, each composite spectrum is further divided by 
the global composite spectrum (blue line in Fig. [T]) and 
normalized at 3000 A. The results are shown in Fig. [2j 

The quasar composite spectra change with Lboi sig¬ 
nificantly, where the higher Lboi quasars have system¬ 
atically harder continua. We calculate the continuum 
slopes (aj /12 and 01 ^ 24 ) for each composite spectrum and 
plot them as a function of Lboi in the sub-panel of the 
figure. The 0^12 and q ;;^24 values are close at high Lboi 
and show similar trends with luminosity. At the low 
luminosity side (log Lboi <47), the change of the con¬ 
tinuum slope with luminosity is very significant, while at 
the high luminosity side (log Lboi > 47), quasar continua 
roughly keep a constant slope (—0.5). 



\ (lOOOA) 


Fig. 1. — Composite quasar spectrum generated from 91,131 
SDSS DR7 quasars (thi s work, blue line) and the result of 
(jVanden Berk et al.l 120011 . red line). The signal-to-noise ratios of 
these two composite quasar spectra are plotted in the sub-panel. 
The dotted lines mark the regions of continuum windows, i.e., 1350- 
1365 A, 2210-2230 A and 4200-4230 A, which are used to calculate 
OLi ,\2 ^^ 1/24 (see the text for details). 

In the UV-Opt band, the quasar continuum can be 
very nicely characterized by a single power law oc 
and is usually f itted from continuum windows. 
In IVanden Berk et alJ ()2001f) . two continuum windows, 



Fig. 2.— Ratios of the composite spectra in different bolometric 
luminosity bins to that of the whole sample. All the composites 
are normalized to the flux at 3000 A. The range of the bolometric 
luminosity bin (in log Tboi) ^-rid the number of contributing spectra 
for each composite are listed in the legend. In the sub-panel, the 
continuum slopes 0 : 1^12 (diamonds) and ai /24 (triangles) are plotted 
against the median bolometric luminosity of each bin. 

The redder spectra of lower luminosity quasars can 
hardly be explained by the accretion disk model, where 
the peak of the big blue bump moves toward lower fre¬ 
quency (longer wavelength) with increasing luminosity 
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liHubenv et alJ 12^01 : iDavis et'aP 120071) . As discussed in 
iDavis et al.l (|200'^b this discrepancy could be alleviated 
by assuming a higher intrinsic dust extinction in low lu¬ 
minosity quasars. Indeed, there is observational evidence 
that the low lu minosity AGNs have higher intrinsic dust 
extinction le.g.. iGaskell et aI1l2004ll . Such a luminosity 
dependent intrinsic dust extinction scena rio is also con¬ 
sistent w ith a receding dust torus model (|SimDSonll2005l : 
iGul 1201311 . However, even if the dust extinction is in¬ 
dependent of the quasar luminosity, quasars with higher 
extinction will also show lower luminosity and redder UV 
spectra. Such a selection effect may also partly explain 
the — Lboi relation we observed. 

3.2. Trend of ai, with redshift 

In this section, we test whether the quasar continuum 
slope shows systematical variance with redshift. As a 
first step, we divide the sample into five redshift bins 
with bin sizes of 1 from z = 0 to 5 and make composites. 
Similar to Fig. ^ we calculate the ratios of the com¬ 
posites to the global one for different redshift bins and 
sh ow the results in Fig. [3l Consistent with the finding 
of iPentericci et all (|2003[ 1. we find little dependence of 
quasar spectra on redshift. 



X (lOOOA) 


Fig. 3.— Ratios of composite spectra from different redshift bins 
to that of the whole sample. 

However, the SDSS quasar sample is roughly a magni¬ 
tude limit sample and thus high redshift quasars are bi¬ 
ased to high luminosity ones. According to the results in 
Section 3.1, if there are no evolutionary effects of quasar 
spectra, high redshift quasars should also have harder 
UV spectra. Therefore, the apparent non-evolutionary 
quasar spectra implies another redshift dependent mech¬ 
anism that could compensate for the luminosity bias. 
To recover this hidden effect, we make composite quasar 
spectra at different redshift bins after Lboi controlled. 

The quasars are grouped into two-dimensional bins of 
Lboi and z. The bin width of Lboi is set to 0.5 dex in 
logarithm while the width of z is set to either 0.1 or 
0.2 according to the number of the quasars available for 
composite. The minimum number of quasars to make a 
composite spectrum is set to 20. The results are shown 
in Fig. m where the spectral slopes are plotted against 


redshift with diamonds {a^i 2 ) and triangles ( 0 ^ 24 )- The 
continuum slopes within the same Lboi bin are connected 
with lines, while different Lboi bins are represented by 
different colors. We find that in any given Lboi bin, 
quasars at higher redshifts have systematically redder 
UV continuum slopes (for both a^i 2 and a,j 24 :)- 
This result clearly identifies the existence of the red¬ 
shift dependence of quasar spectral slopes and implies 
that the previous debates are probably driven by red¬ 
shift selection bias. 



Fig. 4.— ce^i2 (diamonds) and a^24 (triangles) of quasar com¬ 
posites in different Lboi bins (different colors) as a function of red¬ 
shift. Dust effective density riuv is fixed to 2 X Mpc“^ and 

two extinction models: GB07 and SMC are drawn in dotted and 
dash lines, ai, is probed by 01,24 below z = 1.5 and by 01,12 above 
2 = 1.5. 

4. MODELING THE REDSHIFT DEPENDENCE OF THE 
QUASAR CONTINUUM SLOPE WITH COSMIC DUST 
EXTICTION 

As we have mentioned in the Introduction, the cosmic 
dust extinction is a very interesting mechanism that can 
make the quasar continua systematically redder at high 
redshift. Therefore, it is quite natural to introduce a 
cosmic dust extinction model to explain the results we 
have obtained in Fig. SI Other possibilities will be 
discussed in Section 5. 

We start from a toy model, which assumes that dust is 
uniformly distributed along the line of sight. Considering 
the intervening dust has a comoving number density of 
n, the optical depth of the cosmic dust to the photons 
with wavelength Aq emitted from a quasar at redshift zq 

t(Ao) = / na{ ■ Xo)Dh dz (I) 

Jo 1 + ^ E(z) 

where Dh = c/Hq, E{z) = H{z)/Ho, and cr(A) is 
the cross section of dust absorbers at wavelength A 
(iMore et al.ll2009tl . We assume that the properties of the 
absorbers do not evolve with redshift and (t(A) is scaled 
by the cross section of the absorber at 5500A, 

cr(A) = CTv • e((A)/e((5500A) (2) 

where el(A) stands for the extinction law. 

We also assume that the intrinsic continuum of quasars 
in the UV-Opt band can be characterized by a pure 
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power-law form . Quasars are put to their 

redshifts zg, and then, their observational spectral con¬ 
tinuum can be rebuilt and slopes can be measured 
after considering the cosmic dust extinction from z = zq 
to z = 0. To mimic the observation, we also fit the con¬ 
tinuum slopes in different wavelength ranges for quasars 
at different redshifts (q;j /12 for zg > 1.5 and a,z 24 for 
ZQ < 1.5). In this toy model, besides the extinction law, 
we have two free parameters, one is the combined pa¬ 
rameter ncTy , which represents the effective density of the 
cosmic dust particles, and the other is the intrinsic 
quasar continuum slope. 

We have little prior knowledge of the template extinc¬ 
tion law of the cosmic dust. As the 2175 A extinction 
bump is rarely seen in other galaxies, an SMC-type ex¬ 
tinction curve is preferred. In this study, we test the 
featureless SMC extinction law first. We choose sets of 
representative values of a^o a-nd show the predicted 
as a function of zg, which are represented by the dashed 
parallel lines in Fig.Ul Here, na^ = 2 x 10~^h Mpc“^ is 
tuned to match the observed a ^24 — z relation. 

As can be seen, our simple cosmic dust extinction 
model can reproduce the observed 0^,24 — z relation at 
the low redshift range (z < 1.5) very well. At high 
redshift, our model (SMC extinction curve and na^ = 
2 X 10 ~^h Mpc“^) predicts a much steeper ^1,12 — z rela¬ 
tion than that is observed. This discrepancy can be easily 
solved by adopting either a lower na^ at z > 1.5 or a shal¬ 
lower extinction curve at a far-UV wavelength. To keep 
our model simple, we test the later possibility. Indeed, 
there are studies suggesting that t he extinction curves 
in quasars are flat ter in the UV (|Czernv et al.l l2004 
ICaske 1 e t al.ll20 04ll. We take the flat reddening curve 
from ( Gaskell fc Benk^l2007l hereafter, GB07), which 
is quite similar to the Milky Way extinction law, except 
for the lack of a 2175 A bump. We show the GB07 dust 
extinction models with the same ncv = 2 x 10~^h Mpc“^ 
as the dotted lines in Fig. IH As expected, the GB07 ex¬ 
tinction curve matches the observation much better than 
the SMC extinction curve at z > 1.5. 

In order to have a more intuitive impression of the 
cosmic dust extinction, we take the GB07 dust ex¬ 
tinction model and calculate the observer-frame extinc¬ 
tion Ay and reddening E{B — V) as a function of 
source redshift. The Ay and E{B — V) are plotted 
in the upper and lower panels of Fig. [5] respectively. 
In this plot, we also show the cosmic extinction con¬ 
strained from a few other measurements for compar¬ 
ison. The dashed line shows the con s tant co moving 
dust extinction model of iMenard et al.l (|2010all , which 
is constrained from the statistical detection of dust red¬ 
dening around galaxies up to large scales. The other 
obse rvational const r aints a re also extracted from Fig. 
9 of IMenard et al.l (j20I0all , which includes the results 
from measuring th e excess scatter seen in higher red¬ 
shift quasar c olors (iMortsell &: Goobaill2003ll . using the 
Tolman test (|More et al.l l2009l) , com bining constraints 
from luminosity distances and Hfz ) (lAvgoustidis et all 
I2OO9II . and using Mgll absorbers dMenard et al.l 120081 
see Appendix C of iMenard et al.l[20I0al for details). Fur¬ 
ther study about the impact of cosmic dus t on s uper¬ 
nova cosmology can b e found in iGorasanitil (|2006f l and 
IMenard et al. f|2010bll . 
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Fig. 5. — Extinction magnitude Ay as a function of redshift 
using the GB 07 extinction model in Fig. The open circle shows 
the result of iMenard et aO Il2010all . where the dashed line shows 
their extrapolation from a constant comoving dust density model. 
The other observation al constraints are extracted from Fig. 9 of 
IMenard et al.l Il2010all . In the bottom panel, the model-predicted 
dust reddening E{B — V) is plotted against 2 . 

5. DISCUSSION: OTHER POSSIBILITIES 

In this section, we discuss whether there are any other 
explanations to the redshift dependence of the quasar 
continuum slopes other than the cosmic dust. 

5.1. Evolution of intrinsic properties? 

In Section 3.1, we have argued that the internal dust 
extinction might be the cause of the a,y-luminosity trend 
that we observed. Could it also be responsible for the 
— z relation? In general, we know that metal and 
dust are accumulated through various processes in the 
galaxy evolution, e.g., dust formation in Type B su¬ 
perno va and asymptotic giant branch stars (jAsano et al.l 
I2OI4II . Therefore, galaxies at lower redshifts will typi¬ 
cally have more dust. This global trend is opposite to 
the Oj/ — z relation we observed for quasars. 

However, we know that quasar host galaxies are atypi¬ 
cal. Many observations have shown the presence of large 
dust mass in high redshift (z > 4) quasars; though, the 
origin of t he large amount of d ust in such early epochs 
is unclear (|Valiante et al.l[201lL and references therein). 
Moreover, there are also clues that the extinction curves 
of redshift (z > 4) quasars might be different from thei r 
lower redshift counterparts le.g.. iNozawa et ^ l20I5ll . 
Despite the variations, we have not yet found any sta¬ 
tistical evidence that quasars at higher redshifts show 
systematically higher intrinsic extinctions, at least for 
the redshift range we probed (0.7 < z < 3.2). Neverthe¬ 
less, we emphasize that such a possibility still opens. 

Besides the possible evolution of the intrinsic dust, 
could other physical properties (e.g., the black hole accre¬ 
tion process) of quasars evolve? Since the black hole ac¬ 
cretion is a local physical process, we do not expect that 
they show significant redshift evolution effects once their 
basic physical parameters are fixed (e.g., black hole mass 
or accretion rate). In this work, we group quasars into 
Tboi ~ z bins and do not further consider the variation 
of other physical parameters. However, we have tested 
that whether, even after constraining quasars to small 
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black hole mass Mbh and Eddington ratio AEdd/-bboi 
ranges (0.1 dex) in the original 2D Lboi — bins, the 
anti-correlation between and redshift still exists. 

5.2. Host galaxy contamination? 

Quasars are hosted by galaxies, but they outshine them 
in most cases. However, for the low luminosity quasars 
at low redshift, the contamination of starlight might still 
be significant, especially at the near-UV band, i.e., a^ 24 - 
For a log Lboi 45.75 quasar, a bright star-forming host 
galaxy with = — 23 contributes ^ 20% of the flux in 
the 2000 ^ 4000 A range. More than that, the redshift 
range of the near-UV band we probed is from z = 0.7 to 
z = 1.1. This is the time when the cosmic star forma¬ 
tion rate (SFR) starts to decrease. Therefore, it is pos¬ 
sible that there are larger contributions of the starlight 
at higher redshift, so that makes the quasar spectra look 
redder. 

To test this possibility, we generate a series of model 
galaxy spectra at different redshifts according to the cos¬ 
mic st ar formation history given in iMadau fc DickinsonI 
(|2014ll and their luminosities are assumed to be = 
—23. We then combine the galaxy spectra at differ¬ 
ent redshifts with the composite spectrum of log Lboi ~ 
45.25 quasars. Again, we measure the a^ 2 i for the syn¬ 
thetic quasar-l-galaxy spectra at different redshifts. We 
find that, globally, ai, 2 A becomes smaller (redder) when 
the galaxy spectra is combined into the quasar contin¬ 
uum. However, on the other hand, we do not find that 
the combined spectra get redder toward high redshift. 
The reason is that, while going to the higher redshifts 
from z ^ 0.7 to 1.1, because of the higher SFR, not 
only does the contribution of the starlight in the near- 
UV band become larger, but the near-UV band spectrum 
of the galaxy itself also becomes bluer. 

5.3. Quasar color selection bias? 

SDSS quasars are selected using a c omplicated color- 
based algorithm (jRlchards et al.l [20021 1 . Quasar candi¬ 
dates at high and low redshift are selected based on differ¬ 
ent color-color diagrams. Therefore, the color-selection 
criteria might introduce biases to the intrinsic colors of 
quasars at different redshifts. For example, is it possible 
that the redder quasar continua at higher redshifts are 
caused by the missing of red quasars at low redshift or 
vice versa? 

To test this possibility, we take the reddest (the last 
magenta diamond in Fig. 2]) composite spectrum and put 
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it at different redshifts (from z = 3 to 0) and calculate 
their observed colors. We do the same thing for the bluest 
one (the first green triangle). We find that the tracks of 
the colors of these two composites along the redshift still 
assemble in the quasar color-color selection area. From 
this test, we argue that there is no significant selection 
bias in the colors of quasars at different redshifts in SDSS. 

6. CONCLUSION 

In this Letter, we use a large non-BAL radio-quiet 
quasar sample from SDSS DR7 to study the redshift 
and bolometric luminosity dependence of quasar contin¬ 
uum slope by means of the composite spectrum. Con¬ 
sistent with previous studies, we find that the UV-Opt 
spectrum becomes harder toward higher bolometric lu¬ 
minosity, and we suggest that the internal dust reddening 
might be the main reason for this luminosity dependence. 

When the bolometric luminosity is controlled, we find 
that the quasar UV continuum slope is anti-correlated 
with redshift (redder at higher redshift). We propose 
that the redshift dependence of the quasar UV slope 
might be caused by the reddening of the cosmic dust. We 
build a simple cosmic dust extinction model to quantify 
the reddening of the quasars. Based on our model, we 
find that an effective dust density na^ ^ 10~^h Mpc“^ 
at z < 1.5 can reasonably reproduce the a ,^24 — z rela¬ 
tion. The a,^i 2 — z relation at high redshift (z > 1.5) can 
also be easily reproduced by using a shallow extinction 
curve or adopting a smaller effective density of the cos¬ 
mic dust. The cosmic dust extinction constrained from 
our model is globally consistent with the constrains from 
other independent measurements, and has reached to an 
unprecedented high redshift (z 3.2). 
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